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Abstract 

We consider some issues of the representation in the Hamiltonian form of two problems 
of nonholonomic mechanics, namely, the Chaplygin's ball problem and the Veselova prob- 
lem. We show that these systems can be written as generalized Chaplygin systems and can 
be integrated by the method of reducing multiplier. We also indicate the algebraic form of 
the Poisson brackets of these systems (after the time substitution). Generalizations of the 
problems are considered and new realizations of nonholonomic constraints are presented. 
Some nonholonomic systems with an invariant measure and a sufficient number of first 
integrals are indicated, for which the question of the representation in the Hamiltonian 
form is still open, even after the time substitution. 

First, we consider certain general results related to the method of integration of nonholonom- 
ic systems, which S. A. Chaplygin jH] called the method of reducing multiplier. We generalize 
this method so it can apply to a broader class of systems, so-called generalized Chaplygin sys- 
tems. In the remaining sections, we apply these results to find explicitly the Poisson structures 
and the isomorphisms with other integrable Hamiltonian systems. 

1. Generalized Chaplygin systems 

Consider a mechanical system with two degrees of freedom, such that the equations of 
motion can be written in the form: 

dt\dqi) dqi 52 ' dt\dq 2 ) dq 2 ?1 ' (1) 
S = ai(q)qi + a 2 (g)<? 2 + b(q), 

where L is a function of coordinates and velocities; this function will be also referred to as the 
Lagrangian of the system. 

For a special form of S with b(q) = we obtain an ordinary Chaplygin system [H]. 
S. A. Chaplygin showed that the equations of so-called Chaplygin sleigh can be reduced to 
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the form JTJ) (with b(q) = 0); the system JTJ) can be integrated by applying the below-stated 
method of reducing multiplier and the solution of the Hamilton Jacobi equation. In jTj the 
authors show that the Veselova system is a Chaplygin system. The Veselova system describes 
the rotation of a rigid body about a fixed point under a nonintegrable constraint: the projection 
of its angular velocity onto the fixed axis is zero. We show below that the system of equations 
for Chaplygin's ball on a plane [H] can also be reduced to the form (JT} (but with b(q) ^ 0). 

We will call the system (0) a generalized Chaplygin system (it should not be confused with 
that from [22 ED], where a different generalization of Chaplygin systems was offered!). 

Chaplygin showed that in the case b(q) = |], the equations preserve their form under 
time substitutions 

N(q) dt = dr 

if N does not depend on the velocities. Let us show that this also holds for equations (JH). 

dq% 

Denoting the differentiation by q\ = —— we find 

dr 

dL 1 dL dL dL 1 dN A , dL 



N% dq t Ndfi dq t dq { N dq^^'dq, 



where L(q, q 1 ) = L(q, Nq'). 
Substituting into (|XJ) gives 



d ( dL\ dL fFt d ( dL\ dL ,- 



dr \ dq' q J dqi ~ dr \ dq' 2 J dq 

= Aro 1 fdNdL dN dL 
S = NS + 
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(2) 

^ dq 2 dq[ dq 1 dq' 2 ) ' 

It is known for an ordinary Chaplygin system 0| that if there is an invariant measure with 
density depending only on the coordinates, one can choose N(q) for which S = 0; hence, in 
terms of the new time r, the system can be written in the classical Hamiltonian form. Let us 
give a generalization of this result for the case of generalized Chaplygin systems of the form 
(pCJ) 3 under the assumption that the Lagrangian is a quadratic function of the velocities (not 
necessarily homogeneous). 

d 2 L 

Theorem 1. Let det . 7^ 0, and let the system @J admit an invariant measure with 

oqi oqj 

density depending only on the coordinates; then there is a time substitution N(q) dt = dr such 
that 

1) the function S , defined by depends only on the coordinates: S = S(q), 

2) in terms of the new time, the equations of motion can be written in Hamiltonian form: 



where ^ 

qi fc=i 



and the Poisson bracket is given by 

{q i ,p j } = 5 ij , {pi,p 2 } = S(q), {qi,q 2 } = 0. (3) 

Proof 

Let us apply the Legendre transformation to the initial system (JTJ): 



dq, 

here, 



dL 

P i = — , H = Y,Pi<ii-L 



q t = — A = -— + —S P 2 = -— - —S 

dPS 1 dqi dP 2 ' 2 dq 2 dPi ' 

S = ai {q) qi + a 2 {q)q 2 + b{q) = A l {q)P l + A 2 {q)P 2 + B{q). 
Using the Liouville equation for the density of the invariant measure p(q) dP\ dP 2 dq\ dq 2 of the 
system (J3j), we find 

*G£-*«) + *G* + *<«>)-* 

since p depends only on the coordinates, each of the brackets should become zero separately: 

i * - AM = o, I * + a i(9 ) = o. 

P % p <9<?2 

1 <9L 

Now we write equation for S 1 , taking into account the relation — 7—7 = Pf. 

N aq\ 

S = (NAM + Pi + (NA 2 (q) ~^) P i + B{q). 

Thus, if we choose N(q) = p(q), then S = B(q), and the first statement of the theorem is 
proved. 

The second statement can be proved with a straightforward verification of the equations 
and of the Jacobi identity. 

Remark 1. The time substitution, or the reducing multiplier, for multidimensional systems is closely 
related to the invariant measure, but nevertheless, as shown in [2H], can differ from it. 

Note also that reduction to the Hamiltonian form is useful for application of the methods of the per- 
turbation theory, introduction of action-angle variables, analysis of integrability and non-integrability, 
etc. 
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The Hamiltonian systems with the bracket are used for description of systems with a 
generalized potential (for example, motion of charged particles in a magnetic field) or systems 
with gyroscopic forces [25j. In this case, the closed 2-form S(q) dq 1 A dq 2 is called the 2-form 
of gyroscopic forces. Locally, fl can be represented as the exact differential fl = dio, uj = 
Wi(q)dqi + W 2 (q)dq 2 and the equations of motions ((2j) can be written in the form of the 
Lagrange-Euler equations 

d_ ( dL w \ _ dL w = 
dr \ dq\ J dq { 

L W = L + W(q, q'), W(q, q') = W 1 {q)q' l + W 2 (q)q' 2 , 

here the Poisson brackets for the new momenta^ = Pi + Wi(q) and coordinates qi are canonical 
(i.e. {qupj} = Sij). 

If the manifold j\4, on which the coordinates qi, q 2 are defined is compact, then the criterion 
for fl to be exact is 

/n = o. 

M 

Thus, if J fl 7^ (i. e. the 2-form is not exact) then the generalized potential W and the 
corresponding Lagrangian and Hamiltonian functions have singularities (so-called monopoles) 
[23J. In this case, sometimes it is said that the global representation of the equations of motion 
in the (canonical) Hamiltonian form is impossible. 



2. Veselova system 

The Veselova system describes the motion of a rigid body with a fixed point subject to 
nonholonomic constraint of the form (uj, 7) = 0, where uj and 7 are the body's angular velocity 
vector and the unit vector of the space-fixed axis in the frame of reference fixed to the body. 
Thus, for the Veselova constraint, the projection of the angular velocity onto a space-fixed axis 
is zero. This constraint is reciprocal to the Suslov constraint [21], for which the projection of 
the angular velocity onto a body-fixed axis is zero. 

In the moving axes fixed to the body, the equations of motion can be written as follows 
®%: 

9U 

Iuj = Iuj x uj + uj + 7 x — — , 7 = 7 x uj, (5) 

cry 

where fj, is an undetermined multiplier, I = diag (Ji, I 2 , J3) is the tensor of inertia, ^(7) is the 
potential energy. The undetermined multiplier /1 can be found by differentiating the constraint: 
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In the general case the equations (JHJ) admit the integral of energy and the geometric integral 

H=±(Iw,u,) + U('y), 7 2 = 1, (7) 
as well as the invariant measure d 3 u> d 3 7 with density 

Pu) = VhJ" 1 !)- (8) 
When {7 = 0, there is an additional integral 

F = (Iu>, la?) - (la?, 7) 2 = |Ia? x 7 | 2 , (9) 

and hence, the system is integrable according to the Euler-Jacobi theorem 0. 

Remark 2. The integral © is generalized when the Brun potential is added 0E]- Some other integrable 
potentials are given in [Tl IT3]. 

Remark 3. The system ©, @ with the constraint (07,7) = was rediscovered in paper ^H] almost 
ten years after In [E], an explicit integration was performed using sphero-conical coordinates. 

Remark 4. The Veselova system and the nonholonomic systems (considered below) describing the 
rolling motion of bodies belong to the class of so-called LR- and L + i?-systems on Lie groups Q] 
Ej. Several results on the existence of invariant measure for such systems are known. We do not 
consider here these general results, especially useful for multidimensional generalizations. Note also 
that the general methods of reduction of nonholonomic systems were examined in many papers, see 
for example |19j . 

Remark 5. Generalization of the Veselova constraint (a;, 7) = d / was considered in (§]. Using 
Chaplygin's method of integration for a dynamically asymmetric ball with non-zero constant of areas 
[6J, the author presented an explicit integration of the equations. 

It was shown in [l] that the Veselova system is the Chaplygin system JTJ) with b(q) = 
and, therefore, upon the time substitution N dt = dr, it can be written in the Hamiltonian 
form, where the reducing multiplier is N = p~ x . Let us show this explicitly using the local 
coordinates (namely, the Euler angles 6, <p, ip) and then apply the obtained canonical Poisson 
structure of the cotangent bundle of the sphere T*S 2 to construct an algebraic Poisson bracket 
of redundant variables a?, 7. With such an algebraization of the Poisson structure, one can 
naturally establish an isomorphism with the Neumann system describing the dynamics of a 
point on a sphere in a quadratic potential. This isomorphism was straightforwardly established 



in [SI [Zj- Later, we will see that this analogy can be directly extended to the Chaplygin ball 
and the general Clebsch system (which includes the Neumann system as a particular case). 
In terms of the Euler angles, the body's angular velocity uj and the unit vector 7 are given 

by 

U3 — (ip sin 9 sin tp+9 cos tp, ip sin 9 cos tp— 9 sin tp, tp cos 0+<p), 7 = (sin 9 sin tp, sin 9 cos tp, cos 9). 

(10) 

The equation of the constraint is 



/ = (o>, 7) = ip + cos 9tp = 0, 



(11) 



Eliminating the undetermined Lagrange multiplier from the equations of motion we can write 
the equation for 9, tp as a Chaplygin system: 

d (dT 



Jt \M ) ~ ~d9 + ~d9 ~^ ' Jt{lhp~)~ lfy + ~£hp~ 



-9S, 



S 



dT n 



dip 



:i2) 



= sin 2 9 y9(I 2 — h) sin tp cos tp — tp(I\ cos 2 ip + I2 sin 2 tp + I 3 ) J , 

ip=— cos 6ip 

where U is the potential energy of the body in an external field, T = |(u?,Iu>) is the kinetic 
energy without the constraint, while T is the kinetic energy from which ip is eliminated using 
the constraint 



T = Tn 



lf!=— COS t 



-Ii(9cos(p — sin c/? sin 6* cos 9) 2 + 



+ ^^(^ sin ip + tp cos tp sin 9 cos 9) 2 + ^h^ 2 sin 4 (13) 

Remark 6. The representation lfT2|) is obtained upon differentiation under the condition pi[): 
dT_dJh dT_dIh CQS0 ^o. OT _ dlh - sin0 ^£ dT _ 9Tb 

Theorem 2 ([lj). After the time substitution N dt = dr, N = (j,!^^ 1 ^ 2 , the equations of 
motion of the Veselova system take the form of the Euler-Lagrange equations: 

d ( dL\ dL 



d fdL\ dL „_ 
d^[d9'J~~d9~ ' rf^V^ 7 / ' 



(14) 



where L = T — U 
written in the form: 



e=N9', ip=N(p> 



is the Lagrangian function; after time substitution it can be 



1 (V x 7,I(Y x 7)) 
L= 2 (7,1-7) ^ 
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Proof 

The proof is based on a simple computation test: after the time substitution, the right-hand 
side S of ((HI), calculated by virtue of (J2J), should vanish. 

The canonical Hamiltonian form of the equations of motion ((Til) can be obtained using the 
Legendre transformation 

dL dL rr I r 

P8=qqp ^ = <^7' H = Pe e +p vV -L, 

de _dH dip _ dH dp e dH dp v dH ^ 
dr dpg ' dr dp v dr 89 ' dr dip 



Using the canonical variables of (|TH|) and the time substitution (p^Vdet I) _1 d£ = dr, one 
can write the equations of motion of the Veselova system in the Hamiltonian form on the 
(co)algebra of the Poisson brackets e(3): 

m = /U 1/2 u,, r = pz 1 i- 1I2 i, 

dM dH dH dT dH ( 16 ) 

^ = Mx aM + rx af' dr = rx aM' 

F=~(r,ir)(M,M) + tf(r), (17) 

where U (T) = U (p^I 1//2 r) , and, respectively, 

y = r 2 = 1, ( Wj 7) = (m, 7 ) = 0, Pw = (7, i- 1 !) 112 = (r, it)- 1 ' 2 , 

{M 4 , Mj} = e ijk M k , {M i} Tj} = e ijk T k , {T h Tj} = 0. 

Thus, we have a Hamiltonian system with Poisson brackets corresponding to the algebra e(3), 
and the four-dimensional symplectic leaf of the structure corresponding to the real motion is 
given by 7 2 = 1, (M, 7) = (for the classical Euler-Poisson equations a similar situation takes 
place if the constant of areas [23J is zero). Note also that the system ((Hp) , (TP?]) determines a 
certain integrable potential system on a two-dimensional sphere and defines thereby a certain 
geodesic flow. 

The inverse transformation is 

u> = (r, irf 2 r l i 2 M, 7 = (r, ir)- 1 / 2 i 1 / 2 r. 

Therefore, a search for integrable potentials for the Veselova system is now reduced to 
the well-studied problem of search for integrable cases in a Hamiltonian system on e(3) with 
Hamiltonian ijTTjl . So, if U = 0, then the additional integral (pj) can be written as 

f = (im , M)(ir, r) - (im, r) 2 . 



(Note that {H, F} = only on the level (M,j) = 0.) 

It was noted in Ej that for U = the system © is equivalent to the Neumann problem. 
As we can see, such equivalence is not a result of the natural reduction of the Veselova system to 
the Hamiltonian form (fTH| . (|T7j) on e(3). It turns out that the isomorphism with the Neumann 
system is caused by existence of a transformation that does not conserve the Poisson brackets 
but reduces the vector field to the required form on the level surface H = const. 

Indeed, let us consider a Hamiltonian system on e(3) (in the case (M, 7) = 0) defined by 
the Hamiltonian 

H = a l -M 2 (TJY)+^ ((M,IM)(r,ir) - (M,ir) 2 ) . (18) 
It is clear that both terms are the first integrals of the system. The following holds true: 

M 2 (r,ir) 

Proposition 1. On a fixed level — = c and (M, T) = 0, the vector field generated by 

det I 

the Hamiltonian (fT8|) is isomorphic to the vector field of the Kirchhoff equations in the Clebsch 
case with zero value of constant of areas (L, s) = 0; this field can be written in the form 

s = k(as x L + Bs x IL), , 

A; = -v / detI. (19) 

L = k (acs x Is + 3(L x IL — c(det Tjs x I _1 s)) , 

Proof 

Let us change the variables 

l = r 1/2 M, s = (r, ir)~ 1/2 i 1/2 r, 

so that the relations s 2 = T 2 = 1, (M, T) = (s, L) = hold. By virtue of linearity we consider 
the two cases a — 1, (3 — and a — 0, j3 — 1 separately. In the first case the equations of 
motion in terms of the new variables read 

s = -Vdetl (^s x L + (s, L) T~j~i^J^ J > 

L = ~>/drtI 1 { ] L ^\ , s x Is. 
detl(s,l~ 1 s) 

Hence, taking into account (s,L) = 0, j — ' . = Af 2 (r,ir) = cdetl, we get the required 

ys, 1 s) 

result. 

The case a — 0, (3 — 1 can be considered analogously. 

If we consider c as a constant parameter, then the vector field (fT9| is generated on e(3) by 
the following Hamiltonian 

H = ka Qm 2 + |(r, ir)) + k[3 Q(M, IM) - | det i(r, ir)^ . 



If a = 1, (3 = 0, we obtain the Hamiltonian of the Neumann case, while at a = 0, (3 = 1 this 
is the Hamiltonian of the Brun problem. With arbitrary a, (3, this Hamiltonian corresponds 
to the general Clebsch case in the Kirchhoff equations [l2l[2H]. Using the representation JHJ), 
(fTTft . we easily obtain the following theorem for the Veselova system: 

Theorem 3 ([5, 7J). After time substitution, the vector field of the Veselova problem (with 
U = 0) on the fixed level of the integral of energy H = h = const becomes isomorphic to the 
vector field of the Neumann problem. 

3. Chaplygin's ball 

Consider the problem of rolling without sliding of a balanced, dynamically asymmetric ball 
on a horizontal plane in the axisymmetric potential force field (we assume that the geometrical 
center and the center of mass coincide) . We fix a moving frame of reference to the body and 
write equations of motion in the following form pj: 

dU 

M = M x oj + 7 x — — , 7 = 7 x w, 

d ~1 (20) 

M = \uj + Dj x (uj x 7) = Iqu, D = mR 2 , 

where uj is the ball's angular velocity, 7 is the vertical unit vector in the moving frame of 
reference, I = diag (ii, J 2 , 1 3 ) is the ball's tensor of inertia with respect to its center, m and 
R are the ball's mass and radius, and U = U(j) is the potential of the external axisymmetric 
field. The vector M is the ball's angular momentum with respect to the point of contact. We 
present the tensor Iq in the form 

Iq = J - D 1 ® 7, J = I + DE. 

Equations (JHJ) (with an arbitrary potential) admit the integral of energy, the geometric 
integral and the integral of areas: 

H=~(M,u) + U(rf), (7,7) = 1, (M, 7 ) = c = const. (21) 

They also admit the invariant measure indicated by Chaplygin pj, p M d 3 M"rf 3 7, with density 

p M = (detI Q )- 1 / 2 = [detj(l- J D( 7 ,J- 1 7))] _1/2 . (22) 

If there is no external field (U = 0), the system ^ has an additional integral 

F = (M,M), (23) 

hence, it is integrable according to the Euler-Jacobi theorem (HJ. In p], the solution of (|20l) 
was given in terms of hyperelliptic functions. 
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Remark 7. The integral (|2"3|) can be generalized to the cases of the Brun field U(j) = §(7,17) QH 
and gyrostat (2EJ- Other integrable potentials (with the zero constant of areas (M, 7) = 0) can be 
found using the representation of the system on algebra e(3), which is given below. 

Comment. Although Chaplygin developed the general method of reducing multiplier he did 
not apply it to the system ifflty . The paper ^H] indicates possible obstructions to the application 
of this method to the system l|2()jl. On the other hand, already in [2] the following question was 
formulated: if the system I120H can be represented in Hamiltonian form? The problem of Hamiltonization 
of Chaplygin's ball was formulated more strictly by Kozlov [T2| and Duistermaat ^2 El- In [E], the 
authors showed numerically that without a time substitution, the equations of motion of the Chaplygin 
ball are not Hamiltonian because the periods of motion for the orbits lying on the two-dimensional 
invariant resonance tori are not the same. However, after an appropriate time substitution the system 
(|2"U|) becomes Hamiltonian and Poisson brackets are found explicitly (see Ref. |2J). 

Unfortunately, the authors of the review |17j did not succeed in an explicit verifying our result 
(perhaps, due to certain misprints in Here, we will prove the result of (2] using another method 
and then reveal an interesting isomorphism between the Chaplygin ball and the Clebsch case in the 
Kirchhoff equations. Another isomorphism has been indicated in |.31| . 

To describe the ball's rotation let us add to the system (ffllj) the equations for the remaining 
direction cosines: 

a = a x w, f3 = (3 x u>. 
Such a system admits two additional integrals linear in velocities 

(M, a) = const, (M, (3) = const. 

Under such an extension, the integral manifolds remain two-dimensional. Hence, the Chap- 
lygin problem with U = is degenerate or, as it is sometimes called, superintegrable. From 
this viewpoint, Chaplygin's ball is a nonholonomic analog of the Euler-Poinsot top, a well- 
known noncommutative integrable system, and the phase space of this three-degree-of-freedom 
Hamiltonian system is foliated into two-dimensional tori (not three-dimensional according to 
the Liouville theorem). 

It was shown in [2] that, for an arbitrary potential, after the time substitution and change 
of variables 

Pil dt = dr, L = p^M, (24) 
the equations of motion (JHJ) take in the Hamiltonian form: 

dM k , d-f k 

—j— = {H, M k }, ^ = {Hn k } 
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with the nonlinear Poisson bracket 



{Li, Lj} = e ijk (L k - D(L, l)p\JiJjlk) , {L { , 7^} = e ijk ~j kl {ji, lj} = 0, (25) 
and the Hamiltonian is the energy ©, which can be written in the form: 

H=^((l-D(j, J"i 7 )) (L, J-*L) + D (J-*!,, 7 ) 2 ) + 17( 7 ). (26) 



Now we show that the system (|20|) describing the motion of Chaplygin's ball is a generalized 
Chaplygin system (|T}, and the bracket |25| can be obtained using the method of reducing 
multiplier (see theorem [TJ. 

As in the Veselova problem, we use now the local coordinates: the Euler angles 8, <p, and 
the Cartesian coordinates of the ball's center x, y. In the moving frame of reference aligned 
with the ball's principal axes, the angular velocity vector and the normal to the plane are given 
by PJ. 

The equations of the constraints (corresponding to the no slip condition at the point of 
contact) can be written in the form 

f x = x — R8 sin %j} + Rp sin 6 cos '0 = 0, f y = y + R8 cos + Rip sin 8 sin ip = 0. (27) 

The equations of motion with Lagrange multipliers are follows: 



d fOT \ d fdT , , 

A *' 17 - A 2/' 17 



dt \ Ox J dt \ Oy 

it fon\ ot df x df y ,1 fu:i: t) \ on 





l = o, 






<9T 

dp 


op 



(28) 



dt\88 J 08 " x 08 ' ' v 08'' dt\0p J Op " x ^ ' dy? 

where T is the ball's kinetic energy without taking into account the constraints l(2*7jl (obviously, 
this energy does not depend on x, y, and 0): 

T = Kn{x 2 + y 2 ) + 

Eliminating the undetermined multipliers A x and Ay with the help of the first two equations in 
(|28|) and the constraints (f27| we get 

A*^& + A v ^k = -mi? 2 (0 + 00 sin 0), 



dfx Of 

\ x — — + A„— — = —mR 2 (p sin 8 + dp cos — #0) sin 0. 

Op Op 



Hence, the equations of motion for the angles 8 and p do not depend on 0, but only on rp. 
Therefore, is a cyclic variable and can be eliminated using the Routh reduction procedure; 
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after that the equations of motion for 9 and ip can be written as 

dt \de) ~^9~ 9 ' ^ ~ ' (29) 

S = mR 2 sin 9 {^p cos + ip j . 

Here, TZ is the Routh function: 

• dT 

K(9,ip,9,ip) = T - V>— y, 

from which x and y should be eliminated using the equations of constraints, while ip is eliminated 
using the equation for the cyclic integral, 

dT 

— - = (Ii — I 2 )9 sin 9 sin ip cos <p + I 3 cp cos 9+ 
dip 

+ ((ii sin 2 <p + I 2 cos 2 tp) sin 2 9 + I 3 cos 2 9) ip = c = const. (30) 

So, we get the equations of motion in the form of a generalized Chaplygin system (JQ, 
and, since the system has an invariant measure, it is possible to write equations ([211 in the 
Hamiltonian form with the bracket iffil. 

Perform the time substitution of the form 

N(8,tp)dt = dr, (31) 

where N = p^ is the density of the invariant measure (|2*2l . 

According to (J2J), the equations of motion, in terms of the new time, are 

£(^)_^_ >q d ( &1Z \ ^ - /?'<? 

d^V^y ~d~9~~ Lp ' H 1 [d^ ) ~ ~dip ~ ' (32) 

S = c(I 3 + mR 2 )mR 2 N 3 sin 9{I X cos 2 ^ + h sin 2 <y? + mi? 2 ), 
where 9' = ^ = iV"^, ^ = ^ = AT 1 *, 6 = 0, p, p)\ e=m , . = ^,. 



By applying the Legendre transformation to the system ()32|) we arrive at 



Theorem 4. Upon the time substitution p^dt = dr the equations of motion (|29l) for Chaplygin 's 
ball can be written in the Hamiltonian form: 

d9 _dH dp e _8H_-dH dp _ OH d Ptp _dH dH 
dr dpg ' dr 89 dp v ' dr dp^ ' dr dp dpo ' 

with the Poisson bracket of the form 

{0,Pe} = {<P,P v } = l, {p^,Pe} = S(9,p), {9,p} = 0, (33) 
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where 



H = 6'pe + <p' P(p -TZ = 



\vl (iJn ~ 0(7, 17)) + +\-^~ e (hh sin 2 9 + lj 12 cos 2 9 - D( 7 , 1 7 )) + 



PeP v 



rl3(Ii — ^2) cos 6* sin 9 sin </? cos y? ; — — (1% — I 2 ) (I3 + D sin 9) sin <^ cos tp- 



sm9 ' sin 1 

A^r-o ~ N 2 c 2 ~ 

—^MI 21 + D) + —{h + Dsm 2 9){I 2l + D), 
sin 9 sm # 



J12 = h sin 2 <f + I2 cos 2 y?, /21 = A cos 2 <p + I 2 sin 2 </?. 

Expressing the variables L = p^M (|2~3| in terms of the local variables 9, <p, pg, and Pv> , we 
find 

cos # sin y? Ar siny3 cos # cosy? cosy? 

L x = pgcosp-p^ — hciV^— , L 2 = -pe sm ip - — hcA^— , L- i =p H> . 

sm 9 sm 9 sm 9 sm 9 



Using such a transformation one can straightforwardly obtain (J2i 

Let us consider in more detail the integrable case U — with zero constant of areas (M, 7) = 
(L, 7) = 0, because the bracket (OKI) in this case corresponds to the algebra e(3). We write the 
Hamiltonian (jUj) (omitting unessential multipliers) and the additional integral (|23| in terms 
of the variables L, 7 as follows 

H = U 2 ( 7 , B 7 ) - \ [(L, BL)( 7 , B 7 ) - (7, BL) 2 } , 
F = L 2 (7,B7), B = 1 -DJ- 1 = 13 \ 
Using proposition [2 from the previous section, we get the following result: 

Theorem 5. On a fixed level (M,M) = const and (M, 7) = the vector field ([20]) . with 
U = 0, after the time substitution and the change of variables 

s = (7, B 7 )- 1 / 2 B- 1 / 2 7 , L = B-^ 2 L 

is reduced to the vector field of the Clebsch case in the Kirchhoff equations with zero constant 
of areas. 

4. Realization of constraints. Chaplygin's ball with the Veselova 
constraint 

In the paper by A. P. Veselov and L. Ye. Veselova [Zj, as well as in [lZj, the authors consider 
the problem of rolling on a plane of a balanced, dynamically asymmetric ball (Chaplygin's ball) 
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with an additional nonholonomic constraint (the Veselova constraint): 

where E is the unit vector of a space-fixed axis. 

The case E J_ -7, with 7 being normal to the plane of contact, was considered in £[] where 
the system was proved to be integrable according to the Euler-Jacobi theorem. Besides, a 
realization of this constraint by means of absolutely smooth walls was offered in p] (see Fig. 
H} so that the ball moves along the straight line on the plane. 




Fig 1: A ball rolling along a straight line (A. P. Veselov, L. Ye. Veselova). 

In [T2j, the case E || 7 was considered and it was supposed that this system should describe 
the motion of a (massive) rubber ball on a plane (an ordinary Chaplygin's ball, considered above, 
was called in this paper a marble ball). Later, we will show that this system is equivalent to 
the Veselova system (0) (and thereby answer the question of the integrability of a rubber ball's 
motion on a plane, formulated in [TT[). 

Here, we consider the general situation, assuming that E is an arbitrary space-fixed unit 
vector. We start with describing a possible realization of such a system (or constraints), since it 
is obvious that the ball-upon-a-plane realization in this case is impossible. Basing on results of 
[213 UH1 E3J EE] , we offer a more general realization for compositions of Chaplygin constraints, 
Suslov constraints and Veselova constraints, using only perfect rolling (without dissipation of 
energy). 

Let us start with a spherical support [20J. This system describes the motion of a rigid body 
with a fixed point O which is enclosed in a spherical shell; the shell touches an arbitrary number 
of massive dynamically symmetric balls with fixed centers (Fig. [2|) . It is supposed that there is 
no sliding at the points of contact of the balls and the shell. As was shown in [2jQ] this system 
is integrable with an arbitrary number of balls. In the case of a single external ball, we have a 
problem equivalent to the problem of Chaplygin's ball. 
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Fig 2: The spherical support (Yu. N. Fedorov) 




Fig 3: Wagner's realization of the Suslov constraint 
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Fig 4: Realization of the Veselova constraint 

Another system, which can be used for realization of nonholonomic constraints, is called 
a nonholonomic joint (H). The original version of this construction was offered in [TH] as a 
realization of the constraint of the Suslov problem (u>, a) = 0, where a is a body- fixed vector. 
In this case, flat wheels (disks) are attached to the body with the fixed point; these wheels roll 
without sliding upon the interior surface of a fixed spherical shell (see Fig. EJ) - It is supposed 
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Fig 5: The nonholonomic joint 



that a wheel is so sharp that its velocity in the direction perpendicular to its plane is zero. 
Similarly, one can consider the motion of a body with a fixed point when the body is enclosed 
in a spherical shell, which is touched by wheels (disks), whose axes are fixed in space (Fig.BJ. It 
is clear that in the simplest case, we obtain the constraint of the Veselova problem (uj,j) = 0, 
where 7 is the vector lying in the disk's plane (this result was mentioned in [TB]). In [Hj, a 
similar realization of a similar constraint was offered: a frame with wheels (disks) is attached to 
the interior or exterior surface of a spherical shell (Fig.lHJ). Such a construction ensures equality 
of projections of the angular velocities uj s of the spherical shell and u>f of the frame with wheels 
(and, correspondingly, of the bodies attached to them) onto the axis E which is orthogonal to 
the plane containing the disks' axes: 

If the frame with wheels is fixed in space, then we have the Veselova constraint, while the 
space- fixed spherical shell gives the Suslov constraint. 

Now let us consider a combination of a spherical support and a nonholonomic joint, where 
the body with a fixed point is enclosed in a spherical shell which is in contact with a single ball 
and a single disk (see Fig. HJ). 

In the frame of reference aligned with the principal axes of the body, the equations of 
constraints read 

Ru> x 7 + #!U>i x 7 = 0, (u,E) = 0, (34) 

where u> is the body's angular velocity, R is the radius of the spherical shell, wj and Ri are 
the angular velocity and radius of the adjoining ball, 7 is the unit vector of the axis through 
the balls' centers, and E is the normal vector to the plane that contains the ball's center and 
the disk's axis. 
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Fig 6: 



The equations of motion with undetermined multipliers look as follows: 

Iuj = Iuj x u> + R-y x N + fiE + M Q , D 1 uj 1 = D 1 uj 1 x uj + R\i x N, 

( 35 ) 

7 = 7x0;, E = E x uj, 

where I = diag J 2) ^3) is the body's tensor of inertia, D 1 is the scalar tensor of inertia of 
the adjoining ball, N = (N\, N 2 , N 3 ) and fj, are the undetermined multipliers that correspond 
to reaction of the constraints (JMJ), and Mq is the moment of the external forces. Using the 
second equation of (|35|) we find that (wi,*/)" = 0, therefore, 

(wi,7) = -(wi,7) = -(wi,7 x w). 

Using this relation and the second equation in (|35|) . we eliminate 7 x AT from the remaining 
equations. As a result, we obtain 

R 2 

Iuj + Dj x(ux7)=Iwxw + ^E + Mq, D = -=~Dx, 7 = 7x0;, E = E x uj. (36) 

Comparing this with (jSJ and (|Z5j) , we conclude that these equations coincide with the equations 
for the rolling of Chaplygin's ball with the additional Veselova constraint; in this case the 
direction of the vector E can be arbitrary. The undetermined multiplier \i can be found from 
the relation (E,u>)' = 0: 

(Jw x uj + M ,\- l E) 
A*="- { EX X E) ' I Q = J -^®^' 3 = I + D. 

By straightforward calculations one can show that if Mq does not depend on uj, equations 
p6|) have the invariant measure p u d 3 uj d 3 7 with density 

Pu} = ((E,lQ 1 E)detI Q ) 1/2 . (37) 
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There are also obvious geometric integrals 

7 2 = 1, E 2 = 1, (7, E) = const. 
dU dU 

In the potential force field Mq = 7 x — h E x and the energy is also conserved: 

H=±(I Q u,u,) + U('y,E), 

where {7(7, -E 1 ) is the potential energy of the external forces. If there are no external forces 
(U = 0) and (E x 7 ^ 0), then there are two additional integrals: 

F 1 = (K,^x 7 ), F 2 = (K,Ex(Exj)), 

(38) 

K = I Q ^-(I Q o;,^)^; 

hence, the system (Effijl is integrable (according to the Euler-Jacobi theorem). 

In order to prove that the integrals F 1; F 2 exist, let us write the equations of evolution of 
the vector K: 

K K x u. 

Hence, the vector K is fixed in space, and all its projections onto fixed axes are conserved, but 
since (K, E) = 0, only two independent integrals remain. 

Therefore, the system ijHfij) (when E x 7 / 0) is almost identical to the system considered 
in [Zj (where an additional constraint (E,~f) = is imposed). Note that the system ijHfill with 
{7 = has not yet been integrated in terms of quadratures. 

In the special case E = 7 the integrals (|38|) are identically zero, but from the equation of 
the constraint (w,7) = we find that (0^,7)' = (^,7) = and, hence, 

Iqu> = iuj — D(uj, 7)7 = Joj. Iquj = Jus — D(Cj, 7)7 = Jo;, lu: x cj = Juj x u). 

Therefore, in this case after the change I —>■ J the system (|36ll is equivalent to the Veselova 
system (jSJ). Thus, we answer the question of integr ability of rolling of a rubber ball formulated 
in E$j. 



5. Nonholonomic Jacobi problem 

Here we consider another class of nonholonomic systems with an invariant measure and 
first integrals; however, it is still unknown, if the systems' equations of motion can be written 
in the Hamiltonian form. 
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Let a dynamically symmetric ball (I = DqE) of radius R roll without sliding upon a fixed 
surface. In a fixed Cartesian frame of reference, the surface, on which the ball's center moves, 
is given by $(£c) = 0. (It is clear that the surface, on which the ball moves, is equidistant to 
the surface $(aj) = 0, but the equations of motion take a simpler form if the radius-vector of 
the center x is used.) 

The equations of the constraint and the equations of motion with undetermined multipliers 

are 

v — Ru x n = 0, n 



V$(x) 



dx ' 



|V$(a:)r 

niv = iv — — . \iui = RN x n, 
where v = x is the velocity of the center of mass, u? is the ball's angular velocity, m is the 
mass of the ball, N is the constraint reaction, n is the normal to the surface, and U(x) is the 
potential energy of the external field. All the vectors are supposed to be projected onto the 
fixed system. 




Fig 7: Rolling of a ball on a surface (G is the center of mass, Q is the point of contact of the 
ball and the surface) 



After eliminating the undetermined multipliers N, we obtain the equations of motion in 
the form 



m + — £ | uj = m(u), n)h + Rn x — — , x = Roj x n. 
R 2 J ox 

where n is expressed, using the equation of the surface, as follows: 

<9 2 $ 1 ^/»\(»\ d 2 $ 



7* Z-f 



dxjdxi 



<T) 2 



- , ^dxij \dxj J dxjdx 

k \ 3 / \ j / j 

Equations (|3H admit the geometric integral and the integral of energy 



x k . 



Mx) 



0, H = -(D + D)u 2 - ^(u,n) 2 + U(x), D = mR 2 . 



Besides, there is also the invariant measure p u d 2 uj d 3 x with density [2 

Pu = |V$(x)| 



(39) 



(40) 



(41) 
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Is it possible (using, say, the method of reducing multiplier) to represent the equations (jHHj) 
in the Hamiltonian jormP. 

Obviously, in the general case the answer to this question is 'no'. Indeed, writing (|39|) in 
terms of the local coordinates on the surface $(x) = 0, we obtain the system of five equations 
with the integral of energy. If this system is Hamiltonian, then the corresponding Poisson brack- 
et is necessarily degenerate (the rank of the Poisson bracket is always even); consequently, there 
should exist a Casimir function F(x,uj) independent of the energy (|IU|) . But if such a function 
existed, it would have been an integral of equations (fHHl) . which, as numerical experiments 
show, generally does not exist. 

There is an important special case of (|39|) . when the surface is a three-axial ellipsoid §(x) = 
(x, B _1 a;) — 1 = 0, B = diag fo 2 , h) (more precisely, an arbitrary quadric surface). This is the 
so-called nonholonomic Jacobi problem similar to the problem of geodesies on an ellipsoid. 
It is shown in |H that in this case there is an additional integral 

(u> x n. B _1 u; x n) , , 

(n,Bn) v ; 

The integral (fl2f is similar to the Joachimstahl integral in the classical problem of geodesies 

k 1 v-^ c- 

and can be extended to the case when the potential U(x) = —x 2 H — > k,Ci = const, is 

2 2 xf 

i 1 

added. 

If the nonholonomic Jacobi problem can be written in the Hamiltonian form (possibly with 
time substitution)? It is a very complicated question. On one hand, as numerical experiments 
show, the system in this case behaves chaotically and is nonintegrable |H], so we cannot find 
here any obstacles for Hamiltonian representation, typical for integrable systems. Besides, the 
two-dimensional Poincare map, which can be constructed on the level surface of the integrals 
(HOI) and (J12J), has a measure and is symplectic; therefore, this map can be generated into 
the flow of a Hamiltonian system (see, for example, (2^j) ■ On the other hand, the method of 
reducing multiplier and the explicit Poisson structure fitting do not give the required results. 
Note that the possibility of hamiltonization essentially depends on smoothness, analiticity, or 
algebraicity of the sought-for Poisson structure. Here, we do not consider these issues. 

Another nonholonomic problem (somewhat simpler than the previous one, because there ex- 
ists a certain integrable limit problem) concerns the system considered in [21]: it is a problem of 
the spherical suspension. It is supposed in this case that a dynamically asymmetric ball (Chap- 
lygin's ball) rolls on the surface of a fixed sphere. An analysis of the problem's integrability is 
given in [2T1I29] . 
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